correlation.trials=rbind(c(0,0),c(0,.5),c(.3,.7),c(.4,.9))
copula.selections = c("t","gaussian")
marginal.selections = c("norm","std")
sd=1
mean=0
nu=3
parameter.selections = list(list(mean = mean, sd =sd),list(mean=mean,sd=sd,nu = nu))
for (g in 1:length(coupla.selections)){
	fo 
			parameters = list(list(mean = mean, sd =sd),list(mean = mean, sd =sd),list(mean = mean, sd =sd))
		}else{
			parameters = list(list(mean = mean, sd =sd, nu = nu),list(mean = mean, sd =sd, nu = nu),list(mean = mean, sd =sd, nu = nu))
		}

		for (i in 1:dim(correlation.trials)[1])){
			p.x2x1 = correlation.trials[i,1]
			p.yx1  = correlation.trials[i,2]
			p.yx2  = 0
			param=c(p.yx1,p.yx2,p.x2x1)	
		
		}
	}
